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We examine whether charged particles injected by self-annihilating Dark Matter into regions 
undergoing Diffuse Shock Acceleration (DSA) can be accelerated to high energies. We consider 
three astrophysical sites where shock acceleration is supposed to occur, namely the Galactic Centre, 
galaxy clusters and Active Galactic Nuclei (AGN). For the Milky Way, we find that the acceleration 
of cosmic rays injected by dark matter could lead to a bump in the cosmic ray spectrum provided that 
the product of the efficiency of the acceleration mechanism and the concentration of DM particles is 
high enough. Among the various acceleration sources that we consider (namely supernova remnants 
(SNRs), Fermi bubbles and AGN jets), we find that the Fermi bubbles are a potentially more efficient 
accelerator than SNRs. However both could in principle accelerate electrons and protons injected 
by dark matter to very high energies. At the extragalactic level, the acceleration of dark matter 
annihilation products could be responsible for enhanced radio emission from colliding clusters and 
prediction of an increase of the anti-deuteron flux generated near AGNs. 


I. INTRODUCTION 

Cosmic rays are detected up to > 10^® eV energies and 
are composed of a variety of particles, which, depending 
on the energy, could be electrons, positrons, hadrons and 
nuclei [1-6]. The existence of these high-energy parti¬ 
cles requires an astrophysical acceleration mechanism to 
ultra-relativistic energies. One popular process is first- 
order Fermi acceleration, in which charged particles are 
injected in magnetized shock regions. Thanks to scatter¬ 
ing on magnetic field inhomogeneities, particles can re¬ 
peatedly cross the shocks and gain energy. This process, 
also known as Diffuse Shock Acceleration (DSA) [7-10], 
is meant to accelerate particles, such as electrons and 
protons, to high energies at a variety of astrophysical 
galactic and extragalactic sites. 

While DSA has been advocated as a possible accel¬ 
eration mechanism in our own Galaxy, the main astro- 
physical site where the acceleration of cosmic rays can 
take place is not yet established. The most plausible 
source appears to be supernovae remnants (SNRs) and 
is consistent with X- and gamma-ray data for electron 
and hadron acceleration respectively [11-13]. However 
the Fermi bubbles could be another powerful accelerator 
in our galaxy. 

Assuming no other source of cosmic rays other than the 
thermal population, the amount of energy required to ex¬ 
plain the observed spectrum could be a possible issue for 
SNRs: protons would take up to 30% of the SNR shock 
energy while electrons would only take 1% [14]. Another 
possible issue lies in the excess of leptonic cosmic rays 


at high energy. Several experiments such as PAMELA, 
Fermi-LAT, HESS, MAGIC, and AMS-02 [15-23] have 
indeed collected leptonic cosmic rays up to around TeV 
energies. However these observations are difficult to ex¬ 
plain with thermally-injected cosmic rays and the stan¬ 
dard DSA mechanism due to energy losses of cosmic rays 
in the Galaxy. Serious consideration of the details of the 
acceleration mechanism of cosmic rays in the Milky Way 
may therefore be useful. 

Similarly, the radio emission from the so-called ‘Tooth¬ 
brush’ relic in the cluster IRXS J0603.3-1-4214 [24, 25] 
and the non-thermal spectra from several AGNs indicate 
that shock-acceleration is also occurring in extragalac¬ 
tic objects and the observed spectrum might also require 
introduction of a new population of cosmic rays. 

Here, we entertain the idea that dark matter (DM) self- 
annihilations sustain a non-thermal source of cosmic rays 
that get accelerated to very high energies (well above the 
DM mass threshold) thanks to astrophysical shocks. This 
hypothesis is justified by the fact that both galactic and 
extragalactic sites possess high number densities of DM 
in their centres, where astrophysical accelerators also are 
located [26, 27], and could therefore explain cosmic ray 
observations. This also solves a long-standing problem 
with injection into shocks, namely how the injected par¬ 
ticles build in energy from non-relativistic velocities. In 
our case the particles are already injected at relativistic 
energies by dark matter annihilations. 

The combination of both DM injection and DSA gives 
rise to an interesting signature. Dark matter injection 
alone gives an energy spectrum (of electrons, protons or 
other particles) bounded from above by the DM mass. In 
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contrast, shock acceleration can bring injected particles 
to much higher energies, producing a power-law distri¬ 
bution over all energies. Indeed in section II, we show 
that when a DM contribution is added to the ordinary 
cosmic-ray component, the power-law behaviour common 
to DSA is maintained, but with a low-energy cut-off set 
by the DM mass. This cut-off is less prominent if the 
energy loss rate is large. In section III we discuss the 
potential for DM injection towards the Galactic Centre 
to create observable signatures in cosmic ray data and 
radio measurements, and in section IV we discuss extra- 
galactic signals. We conclude in section V. 


II. EMPIRICAL MODEL FOR DARK MATTER 
INJECTION AT SHOCKS 


In this section, we present the model that we im¬ 
plement to describe the shock acceleration of particles 
injected by dark matter. Our mechanism is based on 
Fermi acceleration. First we assume that cosmic rays 
are steadily injected by DM annihilations (see Eq. Al). 
Next we assume that a fraction of these particles, e, ends 
up in the shock region, where they get trapped by the 
magnetic field B. The time-scale over which they are 
expected to stay in the acceleration region is defined by 
Ttrap E/ZeBcVs, where Ze is the electric charge of the 
trapped particle [28] and Vs is the speed of the shock. 
During this time the particles may gain energy by cross¬ 
ing the shock and also lose energy through standard pro¬ 
cesses, namely inverse Compton and synchrotron losses 
for the leptonic components and pion or e^e“ pair pro¬ 
duction for the hadronic part. 

Since the shock occurs within a much smaller volume 
than that characteristic of the injection of dark matter, 
the efficiency factor e can be as small as ~ 10“® [8]. This 
small value could be however compensated by a large 
DM number density near the acceleration site. Besides, 
even such a very small value can give rise to the observed 
electron-proton ratio via DSA [29]. 

In this work, we will assume a delta-function for the 
spectrum of electrons or protons injected by the DM into 
the shock region. This is equivalent to assuming that the 
spatial diffusion can be neglected and that the trapping 
of cosmic rays is much faster than the energy losses. 

The physics inside the shock can be modeled roughly 
using the Fermi mechanism. We give the details of our 
empirical model in appendix A. The evolution of the par¬ 
ticle spectrum (per unit volume), dn/dE, over time t is 
modeled by taking discrete time-steps At. The spectrum 
of particles still in the shock region by the end of the 
time-step At is given by 




FIG. 1: Spectra of electrons which have been injected into a 
region of Diffuse Shock Acceleration by self-annihilating Dark 
Matter. Upper panel: Dependence of the electron spectrum 
within the shock region on the Dark Matter mass. Lower 
panel: Variation of the electron spectrum with the energy 
loss rate, assuming a DM mass of 10 GeV. 
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where the first term represents the particles which remain 
trapped without crossing the shock and the second gives 
the spectrum of the accelerated particles. The term &ioss 
is the total loss rate, and the factors (3 and 6 the probabil¬ 
ity for a particle to return to the shock, and the fractional 
energy gain per shock crossing, respectively. These three 
parameters depend on the nature of the particles and the 
properties of the acceleration region (see App. A). 

After many time-steps, the number of particles inside 
the shock region reaches a steady-state: the number of 
particles that escape and those injected compensate. For 
fixed f3 and 6 the form of the resulting steady-state spec¬ 
trum depends on the injection spectrum /DM(E,m) and 
the energy loss rate 6ioss- Since we took a delta function 
for /dm(A, m), we expect all the annihilation products 
injected by DM to have their energy shifted above the 
DM mass threshold. 
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In figure 1 we show the electron (and positron) spec¬ 
trum injected by DM annihilation after DSA. While the 
initial e"'"e“ spectrum at injection has a cut-off above the 
DM mass, the effect of DSA is to significantly alter the 
spectrum by accelerating the cosmic rays above the DM 
mass threshold. As a result, the final (accelerated) spec¬ 
trum contains a low-energy cut-off fixed by the DM mass 
and, above this cut-off, the spectrum follows a power- 
law (see the upper panel of figure 1), consistent with 
mono-energetic injection scenarios such as considered in 
ref. [30]. 

Our results have been obtained using a simple empiri¬ 
cal simulation. We have assumed that the physics of the 
acceleration mechanism does not depend on the exact lo¬ 
cation of the particle injection and have also neglected 
the spatial dependence of the energy-loss rate. Even so, 
the broad features that we investigate in detail in the 
following, such as the low-energy cut-off around the DM 
mass, should remain in a more robust simulation. 


III. ACCELERATED COSMIC RAYS IN THE 
GALACTIC CENTRE 
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Hence for an electron with energy of 1 GeV the aver¬ 
age life-time is ^ 6 • 10® years. For protons we assume 
that the energy losses are dominated by the production 
of e"'"e“ pairs (from proton interactions with photons) 
for which the energy loss rate is roughly three orders of 
magnitude smaller. Additionally we need to make an as¬ 
sumption on the distribution of supernovae around the 
galactic centre and on the spread of the Mach number A4 
of supernova shocks. For our estimate we assume that all 
supernovae shocks have a Mach number of AI « 10 [34]. 
The assumptions made here are likely to be too simple 
and could lead us to overestimate the flux. However we 
are only interested in showing how this mechanism could 
explain the observed high energy cosmic ray spectrum in 
the galaxy. 


B. Acceleration by SNRs 


We now discuss the case of cosmic rays injected by DM 
annihilations in the Milky Way centre. We will consider 
two acceleration sources, namely SNRs and the Fermi 
bubbles. 


A. Acceleration mechanism 


We will focus on the galactic centre where the DM en¬ 
ergy density is expected to be the largest. The number 
(and properties) of SNRs near the centre are unknown 
so the results highlighted in this section are only for il¬ 
lustrative purposes. However these results allow us to 
determine which of the Fermi bubbles and SNRs are the 
most powerful accelerator of electrons and protons. 

For electrons, we assume that the major energy loss 
mechanism is synchrotron radiation due to the strong 
magnetic field (we assume that the loss rate from inverse 
Compton scattering is « 2.5 • 10“^^GeV“^s“^ due 
to scattering off the Cosmic Microwave Background and 
interstellar field components [31] and thus neglect it with 
respect to synchrotron losses^). As an illustration, we 
take the magnetic field near the galactic centre to be 
B = 50 fiG based on the upper limit derived in [32]. 
Hence assuming that the shock speed is Ug « 3 • 10® m 


^ Such an assumption is only justified if the magnetic field is strong 
enough but large values are justified for the sites we consider. 


To estimate the number of DM-injected products 
which could be accelerated by shocks in SNRs, we as¬ 
sume that the average SNR shock sweeps out a sphere 
of radius 30 pc and lasts for ~ 10 years with a rate of 
supernova explosions ~ 0.01 per century [32]. This gives 
an effective total volume for the supernovae of Vsn 
(47r/3)-(30pc)®-0.1 centuries -0.01 SN/century ~ 100pc®. 
Such a volume has to be compared with that for DM an¬ 
nihilations in the Galactic Centre (GC), VqCj which we 
define as being contained in a sphere of radius 1 kpc cen¬ 
tred on the GC. 

The fraction of cosmic rays injected in the acceleration 
region is estimated to be. 


R 


SN 

DM 


eB 


Vsn 

.Vgc. 





( 4 ) 


where e is the efficiency factor accounting for the pro¬ 
portion of cosmic rays trapped in the acceleration re¬ 
gion [25] and B is the boost factor accounting for an 
increase in the DM energy density (pdm) with respect to 
a Navarro-Frenk-White (NFW) profile [35]. The num¬ 
ber density udm is related to the energy density by 
Pdm = ?^DM wdm- 


C. Acceleration by the Fermi bubbles 

The so-called Fermi Bubbles are another potential ac¬ 
celeration site in the Galactic centre. Although they may 
be related to the central supermassive black hole activ¬ 
ity, their detection seem to indicate a large-scale region 
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FIG. 2: Spectrum of protons injected by Dark Matter added 
to the expected spectrum from standard cosmic ray sources 
into a region of shock acceleration at the Galactic Gentre 
compared with data from KASCADE [36], ATIG-2 [37] and 
GREAM [38]. We assume a DM mass of 100 GeV and an 
injection rate with efficiency times boost factor eB ~ 10^ for 
supernovae or eB ~ 100 for injection into the Fermi bubbles. 
The high-energy cut-off is due to the finite trapping time of 
the particles in the shock region. 



FIG. 3: Gomparison of data from WMAP [41, 42] (the so- 
called ‘WMAP haze’: a radio emission detected towards the 
galactic centre) to synchrotron emission from electrons in¬ 
jected into regions of shock acceleration corresponding to the 
Fermi bubbles by dark matter, under two different assump¬ 
tions for the mass and efficiency times boost factor eB of the 
latter. 


of shock acceleration of hadrons and/or leptons [39]. 
Since, they also coincide with the region where the DM 
energy density is the highest, cosmic rays injected by DM 
self-annihilations could get a boost and contribute to the 
gamma-ray emission from the GC, at high energy. This is 
particularly interesting if DM is relatively light 1 — 10 
GeV) [40] as one might be able to explain both the GeV 
excess observed in the Fermi-LAT 7 -ray data and the 
excess or spectral hardening of cosmic rays reported at 
higher energy. 

Assuming that the Fermi Bubbles represent a ^ 5 kpc 
size region of shock acceleration, then the DM annihi¬ 
lation volume is increased accordingly and the injection 
fraction enhanced, following: 
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where Vfb represents the volume subtended by the Fermi 
bubbles and eB is a boost factor representing our uncer¬ 
tainty on the acceleration mechanism multiplied by the 
injection efficiency e. The maximum energy to which in¬ 
jected electrons can be accelerated is usually estimated as 
Emax ~ eVgBRs, where Rs is the scale-size of the shock. 
In this case taking i?s ^ 5 kpc for the Galactic Centre as 
an estimate of the Fermi bubble scale, see below, we find 
Eraayi ~ 100 TeV. 


D. Expected signatures 

We show in figure 2 a comparison between the pro¬ 
ton flux after acceleration and data from various experi¬ 
ments. Whether DSA occurs in SNRs or near the Fermi 
bubbles makes no qualitative difference. In both cases, 
our results indicate the existence of a spectral feature. 


However this feature is only visible if eB ~ 100 in the 
case of the Fermi bubbles and eB ~ 10^ in the case of 
SNRs. Hence we conclude that both are able to accel¬ 
erate protons injected by DM to high energy but the 
Fermi bubbles are a more powerful particle accelerator 
than SNRs when it comes to DM annihilation products. 

Accelerated annihilation products could also lead to a 
radio signature if the acceleration mechanism is efficient 
enough. We estimate the synchrotron power near either 
SNRs and the Fermi bubbles using the expression [31], 


dW 
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where P{v., E) gives the amount of synchrotron radiation 
per injected electron and dN/dE is the steady-state spec¬ 
trum of injected electrons in the shock region. We com¬ 
pare the results to measurements of the ‘WMAP haze’, 
where an anomalous radio emission was detected towards 
the galactic centre around frequencies of about 30 GHz. 

The WMAP observations are shown in figure 3, along 
with the expected radio signature from DM. While DM 
can explain the intensity of the WMAP haze, we note 
that it does not reproduce the haze’s spectral shape and 
cannot therefore be the sole explanation to this anomaly. 
Besides, whatever the DM mass, the value of eB that is 
needed to make the DM signal visible remains very large. 

Furthermore, we note that DM particles heavier than 
~ 10 GeV could generate a feature in the radio spectrum 
that may be visible at higher frequencies and could there¬ 
fore be of relevance for the Planck experiment [41], which 
also feature a haze with a similar morphology [41, 42]. 
However the dominant dust contribution may prevent 
signal extraction of a possible DM contribution at high 
frequencies. Adding polarization to the predicted tem¬ 
plate might nevertheless help to find further evidence for 
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an ’’anomalous” synchrotron component. For the com¬ 
position of cosmic rays in the galaxy around and above 
1 TeV there will also likely be a sub-dominant composi¬ 
tion change of the cosmic rays if there is a bump in the 
spectrum as in fig. 2 from the (highly boosted) acceler¬ 
ated DM. This might for example manifest itself as a p 
contribution to the AMS-02 flux. 

For the specific case of the Fermi bubbles there is an 
additional interesting point in that the cosmic rays only 
diffuse through the interstellar medium at a fraction of 
the Alfven speed. Indeed the diffusion velocity will be less 
than 100kms“^ even in the hot phase of the interstellar 
medium, and thus it will take more than 10® years for the 
cosmic rays accelerated in the Fermi Bubble outburst to 
diffuse to the solar system. Thus the wave of diffusing 
DM accelerated cosmic rays outlined in this paper will 
reach the Earth later with a significant delay. Indeed the 
accelerated protons and anti-protons propagate outward 
towards us in the Galaxy disk and halo as a wave. When 
this wave interacts with e.g. molecular clouds, and if 
there is an increased energy bump and flux at 1 TeV, then 
diffuse gamma ray emission from the cloud interaction 
with the CRs will also be generated, although this might 
not be easily distinguishable from the ambient flux. 


IV. EXTRA-GALACTIC SIGNALS FROM 
CLUSTERS AND AGN 

We now focus on cosmic rays injected in extragalactic 
sites where DSA is also meant to occur. We will discuss 
two acceleration sites: AGN jets [10], where the dark 
matter number density is high [43-45], and merging of 
clusters of galaxies [25]. 


A. AGNs 

While we do expect electron and proton acceleration in 
these sites, we will focus in this subsection on the accel¬ 
eration of exotic particles near AGNs. The acceleration 
of DM induced cosmic rays that do not have primary as- 
trophysical sources, such as anti-deuterons [46, 47], could 
provide indeed a unique signature of the acceleration of 
DM annihilation products by shock acceleration. 

It is worth bearing in mind however that measurements 
of the anti-deuteron flux at extra-galactic energies 
10^® eV) do not currently exist [46]. Even so we can 
estimate the fluxes by re-scaling the anti-deuteron flux 
expected from the galactic centre see ref. [46]. We 
then need to account for i) the potential increase in DM 
energy density around AGNs (Pagn) compared to the 
galactic centre Pqq, ii) the different volumes associated 
with these objects (Vagn for the AGN versus Vgg for the 
galactic centre) and iii) the longer distance to the source 
(f^AGN for the AGN distance versus dec for the galactic 
centre). These different rescaling factors lead to a flux of 


accelerated cosmic rays near AGN jets of the order 
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bearing in mind that such a simple estimate does not 
account for the efficiency of the AGN shock and other, 
potentially important, factors. 

Taking dac ^ 8 kpc, cIagn ~ 10 Mpe, 

^ 100 GeVern”®, p5,gn ~ 10^® GeVern”®, Vagn ~ 
10 “^^ kpc® and Vgg ~ 1 kpc® we find an extra-galactic 
anti-deuteron flux of < 10“^ m“^s“®sr“®GeV“®, 

based on 4 )q§ ^ 10“®^ m“^s“^sr“®GeV“® for 10 TeV 
Dark Matter. Even though this upper bound is likely 
to be overestimated, we have shown that in principle it 
should be possible to detect extra-galactic anti-deuterons 
at sufficiently high energies from DM injection near 
AGNs. 


B. Merging clusters 

Merging galaxy clusters also show evidence of particle 
shock acceleration. A prominent example is the existence 
of so-called ‘radio relics’, which are extended Mpe sized 
regions of radio synchrotron emission, present towards 
the outer edges of merging clusters. The radio emission 
from these relics is strong. It has a magnitude of order 
~ 10^^^ Watts/Hz at 1.4 GHz [25] and associated spectra 
show a clear power-law behaviour, as would be expected 
from shock acceleration of electrons. Further evidence 
is found from the radio spectral index spatial distribu¬ 
tion [48], which steepens on either side of the radio relic 
in some cases. 

Radio relics are typically observed to have low Mach 
number shocks of Mpe scale near the cluster virial ra¬ 
dius. This poses an injection problem as shocks can¬ 
not easily accelerate thermal electrons to sufficiently high 
Lorentz factors [49] . Another problem for standard ther¬ 
mal plasma injection into DSA is the lack of substantial 
hadronic acceleration signals. This suggests a ratio in 
the number of accelerated electrons to protons that is of 
order O.I, very much higher than the canonical value of 
O.OI observed in galactic cosmic rays or the even lower 
values suggested for SNR acceleration; however there is 
no evidence for this ratio from the Fermi gamma ray con¬ 
straints [8]. A natural resolution might be dark annihi¬ 
lation debris injection which naturally gives comparable 
numbers of relativistic protons and electrons, and meets 
the Fermi constraint. The maximum energy to which in¬ 
jected electrons and protons can be accelerated is of order 
7?max ~ eBRs, where Rg is the scale-size of the shock. 
In this case taking Rg ^ 1 Mpe for the merging cluster 
shocks we find Emax 3 PeV. 

This strengthens the case for a more exotic injection 
mechanism, such as DM self-annihilations. Although 
the number of injected particles is smaller than the am¬ 
bient population, cosmic rays injected with an energy 
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FIG. 4: Radio spectrum due to synchrotron emission from 
electrons injected into a region of diffuse shock acceleration 
by DM particles with a mass of 10 MeV and an annihilation 
cross section of (crw) « 10“^®cm^s“^ (and an efficiency times 
boost factor eB ~ 10®), compared with a power-law from the 
pre-existing ‘thermal’ population. We also show the spread 
of measnred values of radio emission for the various clusters 
quoted in ref. [25]. 


of about ^ 10 MeV (corresponding to a DM mass of 
about 10 MeV, which motivates this particular choice of 
DM mass) can potentially rehabilitate shocks with a low 
Mach number. Since dark matter is expected to exist in 
abundance in colliding clusters, it may then be possible 
to observe an enhanced synchrotron radiation emission 
from DM-injected electrons in these relics. It is worth 
noticing nonetheless that tight bounds on the DM self- 
annihilation cross section (cru) have already been set, us¬ 
ing radio emission from clusters [50, 51] and one cannot 
arbitrarily increase the boost factor or the annihilation 
cross section. 

We base our calculations on the so-called ‘Toothbrush’ 
relic in the cluster IRXS J0603.3-1-4214 [24, 25]. We take 
a magnetic field of B = 10 /iG and assume that electrons 
can originate from Dark Matter self-annihilation over a 
volume of 1 Mpc^. Assuming that the DM is distributed 
around the cluster centre according to a Hernquist pro¬ 
file [24], and choosing a DM mass of 10 MeV and an 
efficiency times boost factor eB ~ 10^, we find a syn¬ 
chrotron power of the order 10^^ Watts/Hz at 1.4 GHz, 
that is of the same order as observations. The resulting 
spectrum is shown in figure 4 where we display both the 
DM and the standard ‘thermal’ electrons contributions. 
Note that the thermal emission is assumed to obey a 
power law [25]. 

While our findings show that a DM signal only becomes 
visible for very large values of the eB factor, this scenario 
could be justified by either a greater DM density than ex¬ 
pected for the Toothbrush cluster, especially towards its 
centre, or a larger acceleration region near the relic itself. 
However, due to its spectral features, radio observations 
at higher frequencies (in particular in the 10® — 10^^ Hz 
range, see figure 4) could constrain such scenarios. Hence 


we expect that experiments such as ALMA [52] could be 
a useful instrument to test the possible shock accelera¬ 
tion of DM annihilation products. Note that DM par¬ 
ticles with a mass of 10 MeV are too light to produce 
protons via self-annihilation, and so the radio emission 
from electrons is enhanced without increasing the size of 
the relativistic proton population. This could explain 
the absence of proton acceleration in such objects, as ob¬ 
served. 


V. CONCLUSION 

Diffuse Shock Acceleration at astrophysical sites is the 
preferred process to accelerate both galactic and extra- 
galactic cosmic rays. The DSA mechanism typically acts 
on the ambient thermal population of electrons and nuclei 
producing particle spectra with a power-law form, as is 
observed in cosmic ray data. Fermi bubbles, galaxy clus¬ 
ters and Active Galactic Nuclei could be favourable sites 
for DSA, as they present shock regions, where the signa¬ 
ture of particle acceleration has been observed. However 
DSA of thermal cosmic rays cannot always easily account 
for very high energy cosmic rays. 

In this work, we have considered the injection of Dark 
Matter annihilation products into shock regions. Using 
a simple empirical model we showed that shock acceler¬ 
ation allows DM injection products to be accelerated to 
energies much higher than the DM mass. The values of 
the efficiency times boost factor eB is central to all of our 
estimates and needs to be very high for the DM annihila¬ 
tion products to contribute significantly to the observed 
cosmic ray spectra. However this may be realistic due 
to a boost in the DM number density near acceleration 
sites. Even in the strongest collisional shock one might 
expect to have in a cluster merger such as the Tooth¬ 
brush, the Mach number is moderate and any dark mat¬ 
ter density enhancement can only be a factor of around 
four. However a much larger density enhancement could 
occur towards the galactic centre or even more plausi¬ 
bly towards supermassive black holes such as those in 
M87 and Gen A, due to the presence of a central mas¬ 
sive black hole [27]. The result may be an enhancement 
of the dark matter density by many orders of magnitude 
due to the presence of a spiked density profile. This may 
therefore account for the large boost factors used in this 
work. The jet in M87 offers an especially attractive en¬ 
vironment for shock-boosting of DM annihilation debris 
[53]. Indeed even without such boosting, a gamma ray 
signal is plausibly detectable from jet-DM scattering in 
the case of Cen A [45] . One can even imagine ultraheavy, 
non-thermally produced, DM particle annihilation prod¬ 
ucts that might be jet-boosted to ultrahigh cosmic ray 
energies in the vicinity of dark matter-enhanced density 
spikes around active SMBH by jet acceleration mecha¬ 
nisms as in [54] . However in all cases our proposed signal 
would also have to compete with the prompt signal from 
DM-annihilation products outside the shock region. 
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More conservatively, however, we find that using sim¬ 
plifying assumptions, protons injected near SNRs in the 
Galactic Centre would lead to an observable signature 
if eB ~ 10 ^ while this factor would be about 100 if the 
origin of the acceleration is the kpc-sized Fermi Bubbles, 
thus showing that Fermi bubbles are an efficient accel¬ 
erator of DM annihilation products. Similar values are 
obtained for the acceleration of electrons but it is worth 
noticing that electrons could in addition induce a poten¬ 
tial contribution to the WMAP haze. 

Injection of cosmic rays by DM annihilations near 
extra-galactic sites where DSA is happening could also 
lead to interesting signatures. The observation of ex¬ 
otic particles at extra-galactic energies (> 10^® eV) with¬ 
out any astrophysical primary counterpart, such as anti- 
deuterons, would be an important signature of DM in¬ 
jection at shocks. In addition electrons injected by DM 
annihilation could give rise to an important contribution 
in radio relics present around merging clusters, for cross- 
sections around {crv) ^ 10“^® cm^s“^ and 10 MeV par¬ 
ticle masses, as can be seen in figure 4. Our modelling 
of the shock acceleration and the DM-injection of elec¬ 
trons and protons is admittedly very simplistic. More 
sophisticated simulations, which include the propagation 
of the DM-annihilation products to the shock regions as 
well as more accurate energy losses and magnetic diffu¬ 
sion processes are required to confirm our values of the 
flux. However we expect the general spectral morphol¬ 
ogy (namely a power-law with a low energy cut-off) to 
remain an interesting spectral feature for future cosmic 
ray searches. 
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Appendix A: Details of the Empirical Model 

In this appendix we list the main steps used to calcu¬ 
late the spectra of electrons or protons after injection by 
Dark Matter into regions of shock acceleration. 

1. Initial spectrum: The cosmic rays produced by 
Dark Matter enter the shock and accumulate ac¬ 
cording to 

ci?7i ciTT- 

— {t + M,E) = —[t,E) + enl,y^{av)fYyM{E,m)M 

(AI) 


where udm is the Dark Matter number den¬ 
sity, (cTf) is the annihilation cross section and 
/DM(A',m) is the normalised spectrum of injected 
electrons or protons which depends on the DM 
mass m. Throughout this work we assume that 
{av) ~ 10 “^® cm^s“^ i.e. the commonly chosen 
‘thermal’ value. We make the assumption that DM 
far from the shock can still inject particles which 
reach the shock region, given a suitable diffusion 
model. Hence throughout this work we integrate 
over the entire DM distribution when calculating 
the injection rate, but correct for the smaller vol¬ 
ume over which shock acceleration takes place. 

In this work we assume that /dm(A, m) oc S{E—m) 
for electrons i.e. that all injection occurs at the DM 
mass. Note that we assume the electron cooling 
time is much shorter than their diffusion time (in 
the case of pairs) but this is not the case for protons 
for which the cooling time by photopion or e“'"e“ 
pair production is probably much longer than the 
diffusion time Indeed as can be seen in figure 1 the 
resulting spectrum depends strongly on the value 
of m. 

2. Energy losses: By the end of the time-step the par¬ 
ticles trapped in the shock, including those injected 
at the start of the time-step, will have lost energy. 
Since we are interested in GeV energies and above, 
we assume that the electrons lose energy through 
Inverse Gompton scattering with diffuse light or 
synchrotron radiation [31], while protons addition¬ 
ally lose energy to e+e“ pair production [55, 56]. 
Hence for each particle with index i the energy 
losses are modelled as, 

E,{t + At) = E,{t) - 6 iossi;"(t)Af (A2) 

where fejoss is the total loss rate in units of 
GeV-Is-b 

3. Shock crossing or escape: The particles only remain 
trapped in the shock region for a finite amount of 
time. At the end of each time-step a certain frac¬ 
tion of particles will leave the down-stream region, 
according to the time-scale for particle trapping 
Ttrap- Hence there are (1 — exp [—Af/Ttrapj) par¬ 
ticles which leave the down-stream region. 

For these particles there are two possible outcomes: 
they can either cross the shock again with a prob¬ 
ability P and return to being trapped, but in the 
up-stream region, or they can escape the shock re¬ 
gion, with probability 1 — /?. For example if a parti¬ 
cle scatters off an irregularity in the magnetic field, 
the angle of deflection determines whether it crosses 
the shock or escapes the shock region entirely [28]. 

For the particles j which cross the shock their en¬ 
ergies are incremented according to, 

Ej{t + At) = Ej{t) + d ■ Ej{t) 


(A3) 


where 6 is the fractional energy gain per crossing. 
Indeed the values of (3 and 5 set the spectral index 
according to a = (ln/3/lni5) — 1. The fractional en¬ 
ergy gain for a non-relativistic shock is In d = U/c 
where U is the shock velocity and c is the speed of 
light. For a shock with a Mach number = 10 we 
have Ind ~ 10“® and so for example to get a = 2 


(close to the expected injection spectral index at as- 
trophysical sources [10, 28]) we need ln/3 ^ —10“® 
i.e. a small energy gain per shock crossing, but a 
high probability that a particle will cross the shock 
multiple times. The spectrum will also steepen 
when energy losses are taken into account. 
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